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Abstract 


A cylindrical one-dimensional model of the Apollo cyrogenic 
oxygen storage tank has been developed to study the effect of strati- 
fication in the tank. Zero gravity is assumed, and only the thermally 
induced motions are considered. The governing equations are derived 
from conservation laws and solved on a digital computer. Realistic 
thermodynamic and transport properties are used. Calculations were 
made for a wide range of conditions. The results show the fluid 
behavior to be dependent on the quantity in the tank or equivalently 
the bulk fluid temperature. For high quantities (low temperatures) 
the tank pressure rises rapidly with heat addition, the heater 
temperature remains low, and significant pressure drop potentials accrue 
(over 100 psia in 2 hours) . For low quantities the tank pressure 
rises more slowly with heat addition and the heater temperature be- 
comes high (as much as 500° K (440° F)). A high degree of stratifica- 
tion resulted for all conditions; however, the stratified region ex- 
tended appreciably into the tank only for the lowest tank quantity. 

The calculations also indicate a significant flux of mass from the 
heater tube for high tank quantities. The results have been compared 
with Apollo 14 flight data. For attitude hold conditions (i.e. no 
spacecraft motions) the calculations are in good agreement with the 
data. 


235 



Introduction 


The removal of the mixing fans from the supercritical oxygen 
storage tanks in the Apollo spacecraft introduces the possibility 
that the fluid in the tank can become stratified. Stratification, 
which results from the temperature gradients surrounding the heating 
element used for pressurization of the tank, can have several possibly 
serious effects. First it can cause erroneous quantity probe 
measurements with the potential for precipitating a mission abort. 
Second, if the low density fluid from the stratified layer enters the 
oxygen supply lines, the potential exists for fuel cell shutdown. 
Finally, if the stratified layer becomes too extensive, a perturbing 
force of even extremely low magnitude can cause depressurization with 
the potential for the existence of two phase fluid in the tank. 

The purpose of this paper is to examine the extent of the strati- 
fication and the potential for depressurization which might exist in 
the supercritical oxygen storage tank. The worst case of no natural 
convection, has been assumed for the stratification calculations. The 
depressurization calculations are based on uniform (completely mixed) 
conditions in the tank. Zero gravity has been assumed throughout 
the analysis. 

The approach which has been taken is to use a simple mathematical 
model which facilitiates a solution but retains the essential features 
of the actual cryogenic storage tank. For this purpose a one-dimen- 
sional cylindrical tank model is used. The governing equations are 
derived from the conservation laws and solved numerically on a digital 
computer. Realistic thermodynamic properties based on information 
from the National Bureau of Standards (ref. 1) and transport proper- 
ties from ref. 2 are used throughout the analysis. 
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Symbols 


A 

C ] 

k 

K 

L 

m 

M 

-+ 

n 


P 

q 

Q 

Qleak 

^input 

^stored 

^net 

r 

t 

T 


v 

V 


a 


area, 

specific heat, j/kg-°K 

thermal conductivity, j/m-sec-°K 

arbitrary constant 

length, m 

mass of fluid, kg 

mass of heater tube, kg 

unit normal to tank surface 

pressure, N/m^ or psia 

heat flux, watts/cm^ 

heat rate, watts 

heat rate due to heat leak at tank wall, watts 

heat rate input to heater tube, watts 

heat stored by heater tube, watts 

net heat rate to fluid (see equation 26 ) , watts 

radius of point in tank, m 

time, sec. 

temperature, °K or °F 

radial velocity of fluid in tank, m/sec. 

3 

volume of tank, m 

velocity vector of fluid in tank, m/sec. 
arbitrary function of time and space 
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V 


del operator 


e 

n 

4 > 

p 

a 

0 


emissivity 

non-dimensional Lagrange space-coordinate 
function defined by equation (27) 
fluid density, kg/m 

— ft 

Stephan-Boltzmann constant, 5.6697 x 10 watts/ 
m 2 -°K 4 

function defined by equation (28) 


Subscripts 

o 

w 

HTR 

COND 

RAD 

b 

(’) 

Dil 

Dt 

A( ) 


initial value 
value at tank wall 
value at heater 
result of conductivity 
result of radiation 
value at a boundary 

dl_l 

dt 


substantive derivative of ( ) 
increment in ( ) 
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Analysis 

Model and Basic Equations 

A simplified tank model which includes the essential features of 
an actual storage tank is illustrated in figure 1. A cylindrical tank 
with a radius of .305 meters (1 foot) is assumed with variations in 
properties in the radial direction only (i.e., the model is one- 
dimensional) . To provide a specified heat input for maintaining tank 
pressure a cylindrical heater tube with a radius of .025 meters (1 
inch) is centered in the tank. A uniform heat leak is assumed to 
exist over the wall of the tank since cyrogenic storage tanks always 
experience such a heat leak. In order to retain the simplicity of 
the model the effect of fluid withdrawal is included by assuming a 
uniform mass flux through the tank wall. 

For the model shown in figure 1 fluid exists inside the heater 
tube as it does in the actual tank. The heater tube in the actual 
tank is provided with numerous holes so that fluid may move in or 
out the tube. In the model two approaches were used to calculate 
the thermally induced flow in and out the heater tube. In one approach 
the tube was considered to be a separate tank. By assuming that the 
pressure in the tube is the same as that in the main tank, the mass 
flux from the heater tube can be calculated directly. In this 
approach no attempt was made to include the effect of the fluid 
leaving or entering the tube on the behavior of the fluid in the tank. 

In the other approach the tube and tank were considered to be one 
system. Since the heater tube mass flux can not be calculated 
directly in this method, the mass of fluid in the tube was computed 
at each time by integrating the fluid density over the volume of the 
tube. The mass flux for the tube was then calculated by differentiating 
the mass of fluid with time. 

In this analysis a one-dimensional model of a supercritical cyro- 
genic oxygen storage tank in zero gravity is considered. For a zero 
gravity environment fluid velocities are generated only by thermal 
expansion and mass removal effects and thus are expected to be small. 
Therefore, viscous effects and inertial forces are neglected. An 
examination of the thermodynamic data for cyrogenic oxygen in reference 
1 shows that for the pressure (700-900 psia) and temperatures (100- 
300°K) of interest the speed of sound varies from approximately 150 
to approximately 900 meters per second. The speed of sound is there- 
fore considerably larger than any expected fluid velocities. Thus 
pressure variations across the tank are negligible with respect to the 
mean pressure. The analysis is therefore based on the equations which 
follow: 


257 



continuity : 


( 1 ) 


+ pv- v = o 

Dt 


momentum: 


Vp = 


o 


( 2 ) 


energy : 


pc dt = _/a£\ DP + V* (kVT) 
Dt \3T/ Dt 


V " qRAD 


(3) 


If pressure gradients are retained then the equations are hyperbolic 
and describe sound waves so that small characteristic times result. 
Neglecting the pressure gradients reduces the system to parabolic 
equations with considerably larger characteristic times and consequent 
savings in computer time. 

To obtain a solution to the basic equations ( 1) — ( 3 ) an equation 
of state (i.e., p = f(p,T)), various thermodynamic properties 

and transport properties (k) are required. For 

ion of state and required thermodynamic properties 
were obtained from the data in reference 1. The thermal conductivities 
were obtained from the data in reference 2. The data from the 
references has been arranged in the form of a table from which desired 
values are found using a table look-up scheme with first order inter- 
polation . 

The radiation transport term has been included in the energy equa- 
tion since the work of reference 3 has indicated the possible importance 
of this mechanism. The absorption coefficient data used in that work 
however, is questionable. At the present time work is being done at 
the National Bureau of Standards, Boulder, Colorado to generate more 
meaningful data. When that data is available the effect of radiation 
energy transport can be evaluated with more confidence. The radiation 
effect has not been included in this analysis. 

Stratified Tank Analysis 

Development o f equation. - As previously mentioned the tank 
characteristics of primary interest in this study are the tank 
pressure fluctuations and thermal stratification resulting from with- 
drawing fluid and cycling a heater to maintain the pressure in a 
specified range. The basic equations will therefore be utilized to 
develop an equation giving the tank pressure as a function of time and 
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an equation giving the temperatures in the tank as a function of space 
and time. With pressure and temperatures the density as a function 
of space and time can then be computed from the equation of state 
(data of ref . 1) . 

In the development of the pressure and temperature equations the 
use of the momentum equation is avoided by using a Lagrangian view- 
point for the substantive derivatives. That is, the substantive 
derivative is written as 


DO 


Dt 


[IP] 

Ldt J particle 


(A) 


so that particles (small volume elements) are followed in the solution. 
In the one-dimensional analysis being used the particles or small 
volume elements can be followed by using the condition that the parti- 
cle mass is constant. Thus 


Po dv o - pdV 

where the subscript "o" denotes values at the initial time and the 
unsubscripted quantities are values at a subsequent time. For one- 
dimensional cylindrical coordinates 


p o r o dr 0 = prdr 


( 5 ) 


where r Q and r represent the initial and subsequent position of a 
particle. The simplicity of using the Lagrangian viewpoint is ob- 
tained by transforming to the initial (r Q ) coordinates for the com- 
putation and transforming back to the actual coordinates (r) to 
analyse the results. 

One possible complication does arise with the use of the Lagrangian 
viewpoint because with mass removal occurring particles are leaving the 
system. At each time the initial coordinate of the particle at the 
tank wall must be determined and particles with larger initial coordi- 
nates must be considered to have left the system. To avoid the problem 
resulting from this continuous change in the coordinate system, a set 
of non-dimensional initial coordinates can be defined. The non- 
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dimensional coordinates are given by 


n 


r 


o 


r 


o ,w 


( 6 ) 


where r D> w ds a f unc tion of time and can be obtained by considering 
the loss of mass from the system in time dt 


r o ,w 


dr o,w 

dt 


m 

2lTLp w r 0 #w 


(7) 


r o ,w 


l 


r o,w dt 


( 8 ) 


To transform partial differential equations from r D ,t coord i na tes 
to Tl,t coordinates the required relations are 


raun 

= i riLii 

U r o J t 

r o,wL 9 D j 

riu.1 

r 3( )i 

L3t Jr Q 

' L3t J n 


t 


r °» w ILJ 

r o,w ^ 


(9) 

( 10 ) 


The pressure equation development is begun by dividing the con- 
tinuity equation (equation 1) by p and integrating over the tank 
volume 


/* I d V 
•'V P Dt 


+ f V- VdV 
V 


O 


(ID 


The Gauss Theorem can be applied to the second integral to obtain 
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/IS Uv 

J p Dt 
V 


/ V*n 


dA = 


o 


Since uniform mass flux over the tank surface has been assumed then 


/ V 




where the convention .is taken that a positive velocity will produce 
a negative mass flux. Thus, 


1 Dp m 

P W dV = ft 


V 


W 


( 12 ) 


If the density is considered to be a function of pressure and tempera- 
ture then 


M _ ( i£\ D£ . dt 

Dt \3Py T Dt ' 3T/ p Dt 


and since Vp = o then D£ _ 3£ which is uniform in the volume. 

Dt 3t 

these conditions in equation (12) it follows that 


Using 


at 


JL -J 

Pw 

1 II 

V Pj 

l 3 tJ 

I dt 

ip Dt (13) 

J 

’ 1 

V P 

UpJ 

l T dV 


is the equation governing the pressure changes. Equation (13) and the 
energy equation (equation (3) without the radiation term) both contain 
3p DT so they are solved simultaneously to give 

at and Dt 
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( 14 ) 


II 

l^ 13 

Jv Ip^d 

3P\ j 

p V* (kVT)] dV 

3t r 1 

Jv P 

(Ue.\ — X— 
'hph 

- [i) P ] 

and 



DT 

iaai T 1 

1 3p\ . 1 

Dt 

\3t| P p2 Cp \ 

i3t) pC p V ' (kVT > 


where equation (14) can be used first to obtain 3p which appears in 
equation (15). 

For one-dimensional cylindrical coordinates 
dV = 27TLrdr 


where L is an arbitrary length chosen to produce a desired tank 
volume and 


V* (kVT) = 1 3/kr 3T \ 
r 3r ' 3r / 


Since it is desired to work in the Lagrangian viewpoint the transfor- 
mations defined by equations (5—10) are used so that equations (14) 
and (15) become 


m f 1 (3p\ 3 f k r 2 D 

|l) dn 
3n / 


3p _ 27TL Pw 01 r 0>w p 0 p2^x/p 3n \ r 0>w r 0 


f Pa^o.wl /d£\ _ T f[3£\ l 2 ) 
% P 2 V P / T P 2 CpL\9T/pl ) 

dn 

(16) 
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(17) 


/3T \ = _ / 3 

£.] T 

9P + 1 

3 , / kr^p 

3T\ 

^3t } particle [c 

T Ip P 2 C p 

3t 2 r 

P o r o r o,w 

\ r o,w r o 

3nJ 


„ ^o,w 3T 


where r|,r and r are defined by equations (6), (7), and (8) 
respectively. To find ’the location of particles equation (5) is 
used to obtain 


r 2T r r 

= f° £ a r 0 dr c 

L 2 -* ri Jr 0f i p 


(18) 


Boundary conditions.- Since equation (17) is a second order partial 
differential equation, two boundary conditions are required. For the 
fluid in the heater tube a symmetry condition is applied at the center 
of the tube 



= 0 


(19) 


and the temperature of the fluid at the heater tube wall is equated 
with the heater temperature. 


T = T 
1 i HTR 


( 20 ) 


For the fluid outside the heater tube, the temperature of the fluid 
at the heater tube wall is equated to the heater temperature 

T = T htr (21) 


and at the tank wall a heat flux (heat leak is imposed) 

<W = - (kg) A„ 


(22) 



For two of the boundary conditions the heater temperature is required. 
To determine its temperature the heater tube is considered to be a 
lumped mass and a heat balance is computed. The elements shown on 
figure 1, are considered to obtain 


^stored Qinput " ^ QcOND ^ QrAD 


M HTR C P’HTR 


9THTR 

3t 


Q INPUT 



inside 
HTR WALL 




outside 
HTR WALL 


A HTR 


e HTR a 



T w ) A HTR 


(23) 


which is the equation used to find the heater temperature. 


Uniform Tank Analysis 


For a uniform tank the temperature and density are constant over 
the tank volume. Equations (14) and (15) can be simplified for this 
case by performing the volume integrals in equation (14) and by 
integrating equation (15) over the tank volume. The results are 


12. = 

3t 


m ~ 

(l£\ 

\ 3T/P 

, Q net (24) 

S(l£' 

t V 3p j 

It 

-i-f 

PCp! 

fie. I 
,\3T J 

| 1 2 / V 

'PJ J 


and 

11 Jl£.\ L_ 12. + _L_ Qnet 

3t \3T/p p2Cp 3t pCp V (25) 


where the Gauss Theorem has been used to write 
/ v V- (kVT)dV =J^(kVT) • nclA 

= E ^COND + ^ ^RAD = ^net 


(26) 
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If one defines 



then (25) and (26) can be written in the simpler form 


■<& = 4- m + M o (29) 

3t V V M aet 


3T „ 0T 3P + ^net 

3t p 8t pC p V (30) 


Numerical Method 


The stratified tank calculations are performed by solving the 
partial differential equations given by (16) and (17) with boundary 
conditions given by equations (19) - (23). Equations (6) - (8) 
and (18) are used to convert the normalized Lagrangian space coor- 
dinate n to the physical coordinates r. For the uniform tank 
calculations the ordinary differential equations given by (29) 
and (30) are solved. The numerical method employed is to replace 
the derivaties in each equation by appropriate finite difference 
approximations and advance the solution in time. The integrals in 
equation (16) are evaluated using the trapezoidal rule. 

Finite difference expressions . - For purposes of illustration 
an arbitrary function of time and space a(t,r) is considered. The 
difference expressions are: 

3a = a(t + At ,r)- a(t.r) 

3t At (31) 
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( 32 ) 


3a _ a(t,r + Ar) - a(t,r - Ar) 
3r 2Ar 


3 2 a _ a(t,r + Ar) - 2a(t,r) 4- a(t,r - Ar ) 

3r2 Ar 2 (33) 


For functions of time only an expression like equation (31) is used. 
At the boundary points the difference expression used for 3a is 

3r 


= a(t,rb) - a(t,rh ± Ar) 

^3ryb +Ar (34) 


where the signs (±) are selected depending on whether the boundary 
is a right hand or left hand boundary. 


If a boundary condition of the form 

_ K 

Wb - 

is specified then (^a\ 


te 

(«l- 


7 1 


is written as 


_ [ ~a(t ,rb) ~ a(t,rb ± Ar) 1 

L +Ar -* 


• f Ar 
Ar 


(35) 


(36) 


Stability . - In a finite difference approach such as that used 
the stability of the computing scheme must be considered. The 
equations of interest as far as stability is concerned are equations 
(16) and (17). A linear stability analysis applied to equation (17) 
yields 


< 

2kr 2 p , 1 / 

l£\ 3R 1 

(37) 

_r 0 2 Po 2c p (Ar)2 p2c p \ 

9t/p at 



i. 
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as a restriction on the time step size. However, the integral nature 
of equation (16) precludes any simple stability analysis of that 
equation. Therefore it was decided to use equation (37) to define 
the time steps for both equations. Because of the uncertainty with 
regard to the stability of equation (16) the results must be checked 
by making calculations for a range of time step sizes and noting 
the changes in the solutions. The results of such a procedure are 
discussed in the section which follows: 


Results and Discussion 

Calculated Results 

Calculations for a uniform tank have been made and the results 
are presented in figure 2. The calculations give the time rate of 
change of pressure as a function of temperature for a range of mass 
flux rates both for heat addition (Q ne t = 125 watts) and for no 
heat addition. The results indicate that for temperatures up to near 
critical values (i.e., near critical temperatures are between 155 
and 165°K) the magnitude of the pressure derivative is relatively 
high and varies significantly with temperature whereas at higher 
temperatures the magnitude is relatively low and the variations are 
only slight. The unusual variations in the pressure derivative in 
the near-critical region reflects the rapid changes in thermodynamic 
properties in this region. These results also show that the mass flux 
rate has a strong influence at temperatures below the near critical 
region but has only a weak influence at higher temperatures. 

Stratified tank calculations for four typical tank operating 
conditions have been made. The tank and heater physical properties 
which were used for all the computations are given in Table I. Input 
values for each of the four cases computed are given in Table II. 
Because of the questionable stability of equation (16) each of the 
cases was calculated with the time step size determined by 1/3 of 
the expression in equation (37) and for values successively smaller 
until the solution became stationary with respect to step size 
changes. The time step sizes determined for stability of cases 1 
and 2 were about an order of magnitude lower than the initially 
computed values whereas the initially computed values were found 
to be adequate for cases 3 and 4. The results for cases 1-4 are 
presented in figures 3-14. The quantities shown are the time histories 
of heater temperature, tank pressure, and heater tube mass flux and 
profiles of temperature and density across the tank for pertinent 
simulated times. The pressure change potentials are indicated in the 
results by the differences between the stratified tank pressure and 
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the associated uniform tank pressure. The degree and extent of strati- 
fication in the tank is assessed by examining the temperature and 
density profiles. 

The time histories in figure 3 show that cases 1 and 2 have the 
same characteristic feature, that is many heater cycles occurred 
during the two-hour simulated time of the calculations. In contrast 
less than one heater cycle was obtained for cases 3 and 4. This 
difference in characteristic behavior is probably due to the fact 
that the initial temperatures of cases 1 and 2 are in the sensitive 
range below the near critical region whereas the initial temperature 
of case 3 is in the near critical region and that for case 4 is 
considerably higher.. It should be noted that the relative insensi- 
tivity of the tank pressure to heat input for the lowest density cases 
(cases 3 and 4) causes the heater temperature to reach quite high 
values (Tjitr > 500°K for case 4) . The very high frequency heater 
cycling for case 1 apparently results when the heater temperature 
becomes stationary ("locked" or "trapped") in the near critical 
temperature region. In this regard it is noted that when the 
heater temperature for case 2 is in the near critical region the 
heater cycles at a relatively high rate but the rate diminishes as 
the heater temperature rises above the near critical region. 

Results for the possibility of a pressure change caused by 
mixing indicate significant pressure drop potentials for cases 1 
and 2 (100 and 120 psia respectively at t = 120 min) , initially a 
pressure drop potential and later a pressure rise potential for 
case 3, and significant pressure rise potentials for case 4 (60 
psia at t = 120 min) . These pressure rise potentials are in dis- 
agreement with analyses of other investigators (see for example 
reference 4) and require further study. 

The heater mass flux rates obtained are presented in figure 4. 

The results presented are for the separate heater tube and tank model 
only. These results show that the heater mass flux rates for cases 
1 and 2 are relatively high compared to the mass flux rate of the 
tank whereas the rates for cases 3 and 4 are relatively low. The 
calculations indicate that fluid leaves the heater tube when the 
heater is on and enters the tube when the heater is off. Thus for 
some tank operating conditions this mechanism may be important as 
far as inducing convective motions of the fluid. 

Heater mass flux rate results considering the tank as one system 
are not available because of excessive computer time required. In 
this approach the heater location moves from node to node because of 
the Lagrangian coordinate system used. As a consequence a fine nodal 
spacing and hence very large computational times were required to 
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avoid numerical difficulties. However, a limited number of calcula- 
tions using this approach have been made. They indicate that the 
heater tube mass fluxes may be as much as a factor of 2 lower than 
those calculated by the separate heater tube approach. The gross 
tank behavior however is not significantly effected. 

Profiles of temperature and density in the tank are presented in 
figure 5 for simulated times which were judged to be pertinent as far 
as illustrating the stratifications obtained. It is seen that the 
degree of stratification is significant in all the cases. For 
example, the density in the area next to the heater is 1/2 to 1/4 
of the density in the outer portion of the tank. However, except 
for case 4, the extent of the stratification is confined to a 
region very close (within about .025 m (1 in.)) to the heater. For 
case 4 the density (and temperature) stratification is significant 
and extends far into the tank. For cases 3 and 4 the increase in 
density in the unstratified portion of the tank over the initial value 
is a result of the work done on that portion of the fluid by the 
expansion of the stratified fluid near the heater. 

Comparison With Flight Data . 

The calculations just described have been compared with Apollo 14 
flight data. The comparisons are of the heater temperature and tank 
pressure and are presented in figures 6-9. Cases 1 and 2 are high quan- 
tity cases (95% and 75% respectively) for times when the spacecraft is 
in an attitude hold condition. Cases 3 and 4 are moderate to low quan- 
tity cases (52% and 12% respectively) for times when the spacecraft is 
spinning for passive thermal control (PTC). Since the calculation 
method presented here considers only thermally induced motions, then it 
should be expected that a better agreement between calculated and flight 
data would be obtained for cases 1 and 2 than for cases 3 and 4. An 

examination of figures 6-9 shows this to be true. 

For cases 1 and 2 the calculations correctly predict the general 

characteristics of the heater cycling and the maximum heater tempera- 
ture. The flight data in figure 6 show the existance of the condition 
previously discussed in the calculated data where the heater tempera- 
ture becomes nearly stationary in the near-critical region and the 
heater begins cycling at a high frequency. The comparison in figure 7 
indicates that for early simulation times the calculations do not agree 
well with the flight data. . For latter simulation times when stratifica- 
tions have accrued the calculations agree well with the data. 

For cases 3 and 4 the calculations are seen to over predict both 
the time for the heater cycle and the maximum heater temperature. 
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In figure 8 note that after the heater is shut-off the pressure does 
not decay as rapidly in the calculations as in the data. Also the 
heater temperature is higher in the calculation than in the data for 
the entire heater cycle. These differences are attributed to the 
fact that the present calculation method neglects the convective 
transport mechanisms resulting from the spacecraft motions. The 
differences in the pressure cycle in figure 9 are not adequately 
understood. However the differences in the heater temperature 
again indicate an effect due to convective energy transport. 


Conclusions 

Calculations have been made of the behavior of supercritical 
oxygen in a storage tank in zero gravity while fluid is withdrawn 
and a heater is cycled to maintain the tank pressure in the range of 
5.9295 x 10^ N/m^ (860 psia) to 6.2116 x 10^ N/m^ (900 psia) . A one- 
dimensional model was used and the thermally induced movemoment of 
the fluid was the only motion considered. Four cases representing a 
wide range of tank operating conditions were considered. Based on 
the results obtained the following general conclusions are made: 

1. For initial temperatures below near-critical, the tank 
pressure oscillates at a relatively high rate and the heater tempera- 
ture remains low whereas for initial temperatures above near critical 
the tank pressure changes slowly and high heater temperatures result. 

2. The results for pressure change potentials have indicated 
that for initial temperatures below near-critical significant pressure 
drop potentials accrue (over 100 psia in 2 hours). 

3. For initial conditions below the near-critical values 
the flux of fluid mass into and out of the heater tube is significant 
and may be important as far as inducing convective motion of the fluid. 

4. A high degree of stratification was obtained for all 
cases considered however the stratified region extended appreciably 
into the tank only for the case with highest initial temperature, 

194°K. 


5. A comparison of the calculations with Apollo 14 flight 
data has indicated that for attitude hold conditions the calculations 
are in good agreement with the flight data but for conditions when 
the spacecraft is spinning the calculations over predict the heater 
temperature and the time for a heater cycle. 
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TABLE I 


Physical Properties Used in Calculations 


Tank Radius = 0.305 m (1.0 ft) 

Tank Length = 0.463 m (1.52 ft) 

Tank Volume = 0.135 (4.75 ft^) 

Heater on Pressure = 5.9295 x 10 6 N/m 2 (860 psia) 

Heater off Pressure = 6.2116 x 10 6 N/m 2 (900 psia) 
Heater Radius = 0.025 m (1. in.) 

Heater Mass = 1.13 kg (2.5 lbm) 

Heater Specific Heat = 460. j/kg-°K (0.11 Btu/lbm-°F) 

Heater Emissivity = 0.3 

Computing Nodes inside Heater Tube = 6 
Computing Nodes in Main Tank = 57 
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TABLE II 


Input Conditions for Stratified Tank Calculations 

Case 1 - Attitude hold, approximately 95% quantity 
T 0 = 110°K (198°R) 

P 0 = 5.9295 x 106 N/m 2 (860 psia) 

Po = 1054.5 kg/m 2 (65.83 lbm/ft 3 ) 

Q input = 123 - 03 watts (419.82 Btu/hr) 

Qleak = 8 * 06 watts ( 27.5 Btu/hr) 

m = -0.00011 kg/sec (0.87 lb/hr) 

£ase 2 - Attitude hold, approximately 75% quantity 
T q = 142°K (255°R) 

P Q = 5.9295 x 10 6 N/m 2 (860 psia) 

P 0 = 837.76 kg/m 3 (52.30 lbm/ft 3 ) 

^ input = 123 - 2<i watts (420.53 Btu/hr) 

Qleak = 0 . 45 watts (22.0 Btu/hr) 

m = -0.000116 kg/sec (0.92 lb/hr) 
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TABLE II (concluded) 


Case 3 - Passive thermal control (PTC) -approximately 

52% quantity 

T q = 157°K (283°R) 

P Q = 5.925 x 10 6 N/m 2 (860 psia) 

P Q = 536.4 kg/m 3 (33.49 lbm/ft 3 ) 

Qinput = m*45 watts (380.30 Btu/hr) 

Qleak = 4.55 watts (15.5 Btu/hr) 

m = -.000155 kg/sec (1.23 lb/hr) 


Case 4 - Passive thermal control, approximately 12% quantity 

T q = 194°K (350°R) 

P 0 = 5.9295 x 10 6 N/m 2 (860 psia) 

p c = 148.73 kg/m 3 (9.28 lbm/ft 3 ) 

Qin p ut = 111*45 watts (380.30 Btu/hr) 

Qleak = ^* 25 watts (14.5 Btu/hr) 

m = -.000107 kg/sec (0.85 lb/hr) 
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Figure 1. - Sketch of tank and heater showing essential features of model. 
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Figure 2. - Time 



Repeating at a high freq. 
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Apollo 14 Data 
Calculations 
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Figure 7.- Comparison of calculations for Case 2 with Apollo 14 flight data. 

Flight data begins at 151 hours GET. 
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Figure 8.- Comparison of calculations for Case 3 with Apollo 14 flight data. 
Flight data begins at 30 hours 10 min. GET. 
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